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We introduce the R-matrix of the gl(N\M) XX-type model; the one of the Hubbard-like 
^ . model is defined by "coupling" two independent XX models. In both cases, we show that the 

R-matrices satisfy the Yang-Baxter equation. We derive the corresponding local Hamiltonian 
in the transfer matrix formalism and we determine its symmetries. 

A perturbative calculation "a la Klein and Seitz" is performed. Some explicit examples are 
worked out. We give a description of the two-particle scattering. 
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1 Introduction 



The Hubbard model was introduced in order to study strongly correlated electrons [1,2] and, since 
then, it has been widely studied, essentially due to its connection with condensed matter physics. It 
has been used to describe the Mott metal-insulator transition [3,4], high T c superconductivity [5,6], 
band magnetism [7] and chemical properties of aromatic molecules [8]. The literature on the 
Hubbard model being rather large, we do not aim at being exhaustive and rather refer to the 
books [9, 10] and references therein. Exact results have been mostly obtained in the case of the 
one-dimensional model, which enters the framework of our study. In particular, the ID model 
eigenvalues have been obtained by means of the coordinate Bethe Ansatz in the celebrated paper 
by Lieb and Wu [11]. 

One of the main motivations for the present study of the Hubbard model and its generalisa- 
tions is the fact that it has recently appeared in the context of N = 4 super Yang-Mills theory. 
Indeed, it was noticed in [12] that the Hubbard model at half-filling, when treated perturbatively 
in the coupling, reproduces the long-ranged integrable spin chain of [13] as an effective theory. It 
thus provides a localisation of the long-ranged spin chain model and gives a potential solution to 
the problem of describing interactions which are longer than the length of the spin chain. The 
Hamiltonian of this chain was conjectured in [13] to be an all-order description of the dilatation 
operator of N = 4 super Yang-Mills in the su(2) subsector. That is, the energies of the spin chain 
are conjectured to be the anomalous dimensions of the gauge theory operators in this subsector. 
In relation to this, an interesting approach to the Hubbard model is given in [14] that leads to the 
evaluation of energies for the antiferromagnetic state and allows one to control the order of the 
limits of large coupling and large length of the operators /large angular momentum. 

There may be the possibility that some integrable extension of the Hubbard model could be 
put in relation to other subsectors of the N = 4 super Yang-Mills theory. Here we will discuss 
a general approach to constructing a number of supersymmetric Hubbard models. Each of these 
models can be treated perturbatively and thus gives rise to an integrable long-ranged spin chain as 
an effective theory. 

Other supersymmetric generalisations of the Hubbard model have been constructed, see e.g. 
[15, 16]. These approaches mainly concern high T c superconductivity models and their relation 
with the t — J model. They essentially use the gl(l\2) or gl(2\2) super algebras, which appear as 
the symmetry algebras of the Hamiltonian of the model. Our approach however is different and is 
based on the QISM framework. It ensures the integrability of the model and allows one to obtain 
local Hubbard-like Hamiltonians for general gl(N\M) super algebras. They can be interpreted in 
terms of 'electrons' after a Jordan- Wigner transformation. 

In this review paper we revisit and slightly extend the results of [17], our goal here being not 
to reproduce the calculations but to focus on the main ideas of our approach. 

The plan of the paper is as follows. Section [2] is devoted to sketch a number of facts for 
the ordinary Hubbard model. In section [3j we define universal XX models. We introduce the 
corresponding Hamiltonians and determine the symmetries of the model. In section^ we summarise 
the construction of the associated Hubbard-like model, in the Shastry and Maassarani approach. 
From the transfer matrix we obtain the Hamiltonian; we also discuss the symmetries. In section [5] 
we perform a second order perturbative computation d la Klein and Seitz [18]. Then, we define the 
Jordan- Wigner transformation, section [6l used in section [7] to give some examples where we write 
explicitly the Hamiltonians in the gl(2\2), gl{4) and gZ(4|4) cases. We finish in section [8] with a 
study of two-particle interactions. 
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2 Hubbard model 



The 1 dimensional Hubbard model introduced by [1, 2] describes hopping electrons on a lattice, 
with an ultralocal repulsive potential that implements a screened Coulomb repulsion, with U > 0. 
The Hamiltonian is given by 

i=i P =t,i i=i 

We will always use periodic boundary conditions. 

In M = 4-SYM theory this model was first observed in [12] , where a magnetic flux <j> of Aharonov- 
Bohm type was included. In that paper, the potential term was written in a slightly different but 
equivalent form. The relation between couplings was identified and the system was taken at half- 
filling; for our needs we just observe that the ratio t/U corresponds to the coupling g 

t Q , 
U = W <2 ' 2) 

We observe that the Hamiltonian is Hermitian if (j) E R. In the following we will work with this 
flux equal to zero. 

The underling algebraic structure leads us to superalgebras: on each site i the fermionic struc- 
ture 

{ C P ,;> c l'j} = 6 P, P 'kj { c P ,n c pij} = H,n c l'j} = ( 2 - 3 ) 

is a realisation of the super-Lie algebra gl(l\l) ffi <?Z(1|1). The full model algebra is obtained by L- 
times the tensor product of the one site structure. We can easily represent the fermionic structure 
by a graded tensor product of Pauli matrices, written here with the standard notation for basis 
matrices E a p to emphasise the grading: 

Ei 2 - P ,i = c p i , E 2 i-p,i = cl ti , E 22 ;p,i = n p i = S p i c p i , £ll ;p ,i = 1 - n p i = c p ^ p i (2.4) 

*>-(!;)■ *.-(;!). Mo? 

When it occurs, the second pair of labels p, i indicates the spin polarisation p and the site i. 
The matrices -E12 , E21 are taken of fermionic character (they satisfy anticommutation relations 
whatever their spin and space labels are) and En , E 22 are taken of bosonic character (they always 
enter commutation relations whatever their spin and space labels are). The relation (12. 4D is a graded 
Jordan- Wigner transformation^ and respects periodic boundary conditional). We now rewrite the 
Hamiltonian in the spin chain language 

L L 

H = — t (^21;p,i Ei 2 - p>i+ i + E 2 i ;p>i+ i Ei 2;Pi i) + U (-E'H;T J i — ^22;T,i) ~ ^22;|,j) 

i=l p=1,l i=l 



1 Th.e ordinary Jordan- Wigner transformation is cj i = i TT "!.* ^ or ^ ne U P polarisation; an additional term 

k>i 

occurs for the down polarisation. 

2 The standard one violates periodicity. 
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and we split it into the sum of the two polarisations 



L 



H - H XX + H XX + U ^2 { E ll;1,i ~ ^22;T,i) { E ll;l,i ~ E 22;l,i) 



(2.5) 



L 



XX — 



t {E21;p,i E 



12;p,j+l + ^21;p,i+l E\2;p,i) ■ 



i=l 



Taking one polarisation of the kinetic term we easily see that 



E21;p,i E\2-p,i+X + ^21;/j,j+l E\2;p,i 




(2.6) 



1 

1 




namely we see the appearance of a (graded) XX spin chain Hamiltoniaro (or better two XX spin 
chains, one for each polarisation) within the Hubbard model. 

It turns out that the breaking of (|2.5p into the Hamiltonian of two XX models plus a potential 
will allow us to generalise this model to higher algebraic structures by maintaining its main property: 
integrabilit}0. 

A first hint of integrability of the Hubbard model came from the coordinate Bethe Ansatz 
solution obtained by Lieb and Wu [11] but a full understanding of it by the existence of an infinite 
set of commuting charges came much later. A complete set of eigenstates was constructed in [19] 
using the SO (4) symmetry of the ID Hubbard Hamiltonian. Within the framework of the quantum 
inverse scattering method, an R-matrix was first constructed by Shastry [20, 21] and Olmedilla 
et al. [22], by coupling (decorated) R-matrices of two independent XX models, through a term 
depending on the coupling constant U of the Hubbard potential. The proof of the Yang-Baxter 
relation for the R-matrix was given by Shiroishi and Wadati [23]. With a standard construction, 
a transfer matrix can be constructed by taking the trace of a tensor product of R-matrices. The 
Yang-Baxter equation guarantees that the transfer matrix is the generating functional of an infinite 
set of commuting charges. One of these charges is the Hamiltonian (|2.5p itself. 

The construction of the R-matrix was then generalised in the gl(N) case by Maassarani et al., 
first for the XX model [24] and then for the gl(N) Hubbard model [25,26]. Within the QISM 
framework, the eigenvalues of the transfer matrix of the Hubbard model were found using the 
algebraic Bethe Ansatz together with certain analytic properties in [27-29]. 

3 Universal XX models 

We generalize the construction given in [17,24,29] to the case of an arbitrary representation space 
V, possibly infinite dimensional. We will use the standard auxiliary space notation, i.e. to any 
operator A 6 End(V), we associate the operator A\ = A ® I and A2 = I <8> A in End(V) (8) End(V). 
More generally, when considering expressions in End(V)® fc , Aj, j = 1,... ,k will act trivially in all 
spaces End(V), but the j th one. 

To deal with superalgebras, we will also need a Z2 grading [.] on V, such that [v] = will be 
associated to bosonic states v £ V and [v] = 1 to fermionic ones. 

3 At this point it should be clear that the difference between graded and non graded cases appears when boundary 
effects are observed; the thermodynamic limit usually ignores such terms, being sensitive to bulk contributions only. 
4 The flux 4> does not affect integrability properties. 
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We will also assume the existence of a (super-)trace operator, denned on a subset of End(V) 
and obeying cyclicity. When V is finite dimensional, dim(V) = K, End(V) is a matrix algebra or 
super-algebra so that the trace operator is the usual trace or supertrace of K x K matrices. When 
V is infinite dimensional, the definition of a trace operator is more delicate and we will just assume 
that it exists and is cyclic, for the operators we use. 

The construction of a universal XX model is mainly based on general properties of a given 
projector and a permutation. Our main projectors are chosen in End(V) as being 

tt:V^W , n = I-TT:V^W with V = W © VV (3.1) 

In the tensor product of two vector spaces we take the (possibly graded) permutation 

|V0V^V0V 

For example, in the superalgebra gl(N\M) a possible choice is 

7T = E jj J 7T = I - 7T = E NN + E N+M ,N+M (3.3) 
jjtN,N+M 

3.1 R-matrix 

From the previous operators, one can construct an R-matrix acting on V (8> V 

#12 (A) = Si 2 Pi 2 + S 12 sinA+(I©I-S 12 )Pi2 cos A (3.4) 
where £12 is built on the projection operators: 

S12 = 7Tl 7f 2 + 7Tl 7T 2 (3.5) 

It is easy to show that Si 2 is also a projector in V (8> V: (Si2) 2 = Si2- 
Let us introduce the operator C: 

C = 7r-n. (3.6) 
It obeys C 2 = I and is related to the R-matrix through the equalities 

£12 = ^(1-^2) and I®I-Ei2 = i(l + CiC 2 ) (3.7) 

In [17] we gave proof of a number of useful properties of the R-matrix. Essentially the same proofs 
work also for the slightly more general formulation given here. The main properties are unitarity, 
regularity, and Yang-Baxter equation (YBE), that guarantees us that we have an integrable model 

-Ri2(Ai2) -Ri3(Ai3) #23^23) = #23^23) #13(^13) #12(^12) 

where Ay = Aj — Xj. (3-8) 

3.2 Monodromy and transfer matrix 

With a very standard construction, from the R-matrix one constructs the (L sites) monodromy 
matrix 

Aki..x>(A) = 22oi(A) ^o 2 (A) • • • R 0L {\) (3.9) 
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where we tensor product one R-matrix for each site of the theory. It obeys the relation 



#oo'(A - M) A<i...l>(A)A)'<i....l>(m) = A)'<i...l>(a) A)<i...l>(A) Rqo'{\ - n) . (3.10) 

where and 0' are two copies of the auxiliary space. This relation allows us to construct an (L 
sites) integrable XX spin chain through the transfer matrix 

h... L (X) = str Aki...l>(A) = str (#oi(A) i? 02 (A) • • • RqlW) ■ (3.11) 

where, if V has infinite dimension, we assume the existence of the supertrace for the previous 
operator. Indeed, the relation (|3.1U|) implies that the transfer matrices for different values of the 
spectral parameter commute 

[ti... L (*),h...M]=0. (3.12) 

Here the cyclicity of the supertrace has been used. 

Since the R-matrix is regular (namely it is a permutation in A = 0), logarithmic derivatives in 
A = give local operators. We choose the first one as XX-Hamiltonian 

H = t^Uor'^iO) (3.13) 
L 

= ^2 H j,j+i witn H j,j+i = p j,j+i 

where we have used periodic boundary conditions, i.e. identified the site L + 1 with the first site. 
After (|3.12[) . we see that any expansion of the transfer matrix in the spectral parameters A, [i 
generates a set of commuting operators. In particular they commute with the Hamiltonian (j3.13j) . 
so are conserved charges. This formally proves that the system is integrable. 

Explicitly, the two sites Hamiltonian corresponding to the example (|3.3p reads 



H jd+1 = Yl J2 {(-l) lj] Ei j ®E ji + (-l)WE ji ®E ij y (3.14) 

i^N,N+M j=N,N+M 

3.3 Symmetries of the universal XX models 



The choice of the fundamental projectors in (|3.ip directly fixes the symmetries of the model. 

One easily shows that an operator M £ End(W) ©End(W) commutes with the projectors (13. 1|) ; 
then it commutes with the R-matrix in the following sense 

(Mi + M 2 ) #i 2 (A) = i?i 2 (A) (Mi + M 2 ) . (3.15) 

Commutation does not hold if the operator mixes the two subspaces. 

As a consequence of (|3.15|) . the transfer matrix also has a symmetry (super)algebra 

S = End(W) © End(W) (3.16) 

with generators given by 

M<i... L > = Mi + M 2 + ... + M L . (3.17) 

The same is true for any Hamiltonian H built from the transfer matrix so (13.171) commute with 
the HamiltoniarJl. 

5 In principle, this construction cannot exclude the existence of operators that commute with the Hamiltonian 
but not with the R-matrix. In that case, these additional symmetries would have the strange feature of not being 
symmetries of at least one conserved charge (by reconstructing the R-matrix from an expansion). 
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We can reverse this construction: we require a symmetry algebra S from which we construct 
the subspaces W and W. This uniquely fixes the fundamental projector tt that immediately leads 
to obtain the XX model possessing S as symmetry. 

The example f|3. 3|) admits S = gl(N — 1\M — 1) ® gl(l\l) as symmetry superalgebra whose 
generators M have the form 

E jk , j,k^N,N + M for gl(N - 1|M - 1) 

E jk , j,k = N,N + M for gl(l\l). ^ °> 



4 Universal Hubbard models 

Starting with universal XX models, one can build universal Hubbard models, in the same way it 
has been done for usual and super Hubbard models [10, 17]. The logic will be to start from two 
possibly different universal XX models of section [3] and "glue" them with the generalisation of the 
construction given in section [2j 



4.1 R- matrix 

We start with the i?-matrices of two universal XX models, R^iX) and i? 12 (A), living in two different 
sets of spaces that we label by f and j. Let us stress that the two XX models can be based on two 
different (graded) vector spaces and V^, with two different projectors ir^ and ir^. 

The Hubbard model is constructed from the coupling of these two XX models. Its i?-matrix 
has two spectral parameters Ai , A2 and reads: 

fl 12 (Ai,A 2 ) = 4 J (Ai 2 )4,(A 12 ) + ^^4 ^Hh[ 2 )Rl 2 (\[ 2 )CiR{ 2 (\' 12 )C± (4.1) 

sm(A 12 J 

where A12 = Ai — A2 and A' 12 = Ai + A2- Moreover, h' 12 = h(Xi) + fo(A 2 ) and the choice of the 
function h(X) is fixed by the proof of the Yang-Baxter equation. Indeed, when the function h(\) is 
given by sinh(2/i) = U sin(2A) for some free parameter U, the R-matrix (|4.ip obeys YBE: 

-Rl2(Ai, A 2 ) i?i3(Ai, A3) i?23(A2, A3) = i?23(A2, A3) i?l3(Ai, A3) i?l2(Al, A2) ■ (4.2) 

As remarked in [17] the proof relies only on some intermediate properties that are not affected by 
the choice of the fundamental projectors (|3.1|) . The proof follows the steps of the original proof 
by Shiroishi [30], in the same way it has been done for algebras in [10]. Moreover, it was already 
noticed in [10] that one can couple two XX models based on different gl(M) algebras: this naturally 
extends to general (graded) vector spaces V. 

The given R-matrix is regular but non symmetric. It satisfies unitarity (we correct here an 
inconsequential typo that occurred in eq. 3.4 of [17]) in the form 

#i2(Ai,A 2 )i?2i(A2, Ai) = fcos 4 (A 12 ) - (^tM tanh(/ l , 12 ) cos 2 ^)) 2 ") Ii ®I 2 (4.3) 

V v sm(A 12 ) J J 

where = I T (g> I 1 . 



4.2 Monodromy and transfer matrix 



We use the construction given in section 
the 'reduced' monodromy matrix 



to obtain the Hamiltonian of the system, starting with 



£o<i...l>(A) = R i(X,fi) . . . R QL (X,fi) 



jU=0 



(4.4) 
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Any other choice for fj, is possible but, at least in view of obtaining a local Hamiltonian, they do 
not give new information. Provided the supertrace exists, the transfer matrix is given by 

ti...L>(A) = str £o<i...L>(A) 

Then, one gets 

[H, t(A)] = , VA , for H = H{0) = t(O)^ 1 1/(0) (4.5) 
The 'reduced' R-matrices that enter in (j4.4H take a particularly simple factorised form 

i? 12 (A,0)= J Rl 2 (A) J R| 2 (A)/ 1 n (/ l ) (4.6) 

where 

'2 



if* (h) =1(8)1 + tanh( J) (4.7) 



and we arrive at a Hubbard-like Hamiltonian 

L L 

H = E H ^ = E p L + i + + u c l c i] ( 4 - 8 ) 

where we have used periodic boundary conditions. 
4.3 Symmetries 

The transfer matrix of generalized Hubbard models admits as symmetry (super)algebra the direct 
sum of the symmetry algebras of the XX components 

S = End(W T ) © End(W T ) © End{W l ) © End^) . (4.9) 

To prove this symmetry, it is useful to remark that (j3. 15[) can be now specialised to the cases up 

and down. Moreover, the up R-matrix commutes with the down generators and viceversa. We also 
check that 

MC" 7 = C CT M, o-=T4 (4.10) 

where 

M = M T + M l and M CT € End(W a ) © End(VV <T ) . (4.11) 

Thus, one gets 

[i?i 2 (A, 0) , M] + M|] = = [i?i 2 (A, 0) , M} + M£] (4.12) 

that can be easily extended to hold for the monodromy and transfer matrices and for the Hamil- 
tonian; the generators of the symmetry have the form 

L L 

M T = J]M] and M^ = ^Mj (4.13) 

j=i j=i 

The ordinary Hubbard case and all the cases where V a is two dimensional are special because, in 
addition to the list of generators contained in (|4.9p . there are new generators given by 



V ± = of®of, W ± = of®of. (4.14) 

To be precise, commutes with the Hamiltonian if L is even while commutes in all cases. 
These additional generators do not commute with H if dim(V CT ) > 2; they are responsible for the 
SU (2) x SU(2) symmetry of the even Hubbard model. 
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5 Perturbative expansion of the Hubbard-like Hamiltonian 

We expand the Hamiltonian (|4.8|) in the inverse coupling y ; according to (|2.2p , this corresponds to 
the small coupling expansion of the gauge theory. Indeed, precisely that expansion has been used 
in [12] to match the SU(2) dilatation operator with the effective Hamiltonian of the Hubbard model. 
The system was taken at half- filling to guarantee the required spin chain behaviour. With the form 
of the potential used in (|2.ip the half-filled condition is enforced by the U — ► oo requirement itself. 

We take the set of all Hamiltonian eigenstates whose leading energy term is —LU, for large 
positive U. These states are selected by the following projector 



that projects on the subspace where, on each site, one and only one among vrj- , 7rj has nonzero 



action. 

We follow the method introduced by Klein and Seitz [18] to obtain an effective Hamiltonian for 
the corrections to the leading energy —LU: 



For L > 2 the second order effective Hamiltonian is 

H S = E H H !j+ i = 2 E( 1+ p h + A^) W ~A -J+i -U+M ( 5 - 3 ) 



3 3 

For the ordinary Hubbard model this expression can be given in terms of Pauli matrices 

L 

r(2) _ 

i=l 



H S = X> - wi+i) ( 5 - 4 ) 



where the fermionic oscillators of (|2.ip have disappeared and only spin degrees of freedom are left 

The structure of the two-sites Hamiltonian H c ^ ii+l can be obtained explicitly. In matricial 
form, it has diagonal block structure, with blocks given by one of the two matrices 

B-={\ ~ l ) or B + =(] JV (5.5) 



1 1 J ^ \1 \ f 

all other entries being zero. The number of appearances of each block depends on the actual model 
under examination. 



6 Jordan- Wigner transformation 

Let us consider p sets of fermionic oscillators cf\cf^ (i = 1, . . . , L and q = 1, . . . ,p) that satisfy 
the usual anticommutation relations 

{C^,cp} = dy V tf,^} = {cP,cP} = (6.1) 



One defines the following matrix (where the usual number operator) 
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The entries xj?a of this matrix have a natural gradation given by [a] + [/?] where [1] = 1 and 
[2] = 0. 

In the gl(2 p 1 \2 P 1 ) case, one defines at each site i the generators 



a-l 



a=2 6=1 

It is easy to verify the following properties: 

( X -i;ai...a p ,a' 1 ...a' p ) X i;a' 1 ...a' p ,ai...a p (6-4) 
- X t;ai...a,,,a' 1 ...a£ X i;/3 1 .../3 p ,j3' 1 ...p p = ^a{Pi ■ ■ ■ ^a' p /3 p X i;ai...a p ,/3' 1 .../3' p ( 6 - 5 ) 
^ ] -Xi;ai...atp,ai...atp = 1 (6-6) 

-Xi;ai...a p ,a' 1 ...c^ X r,f3- L ...f3 p ,l3' 1 ...p> p = (- 1 ) 9 X j;/3 1 ...[3 p ,p' 1 ...l3> p X i;c n ...a p ,a' 1 ...a' p (« / j) (6-7) 

where 5 = (E (M + Kl)) (E fl^l + [$] 

a=l 6=1 

This means that the operators X i- ai ... ap ,a' 1 ...a' p built out of fermionic oscillators are actually a real- 
isation of the gl(2 p ~ 1 \2 p ~ 1 ) superalgebra. A generic case gl(N\M) can be understood as contained 
in the smallest superalgebra for which N,M < 2P . The unwanted states can be consistently 
projected out. 



7 Examples 

It is possible to construct examples of both XX and Hubbard-like Hamiltonians. Clearly, the XX 
ones are "quasi- free models" because they do not contain external potentials and, if written with 
fermionic oscillators, they only contain hopping terms. In spite of this, they show curious "screening 
effects" namely particles that are allowed to move only if particles of other types are present (or 
absent, depending on the case). We will concentrate on universal Hubbard model examples. 

The first example to cite is, of course, the original Hubbard model of section EJ that is described 
in this formalism as (?Z(1|1) ® gZ(l|l) with the choice tt^ = 7r^ = En, tt^ = tt * = £22. 



7.1 gl{2\2) © gl{2\2) Hubbard Hamiltonian 

This is a more complete example of the models under examination. It precisely implements two 
copies (up and down) of the example (|3.3I) with N = M = 2. The kinetic term of the Hamiltonian 
has a factorised form 

L 

Huub = E { E ( C U C a,i+l + C l,i+l C a,i) { C a/a,i+l + C 'l,i+A,i + 1 ~ K,i ~ 
*=1 <r=U 

+ U(l-2n hi )(l-2n Ui ) } (7.1) 

where the factor 

K,i,i+i = i c aA, i+ i + c *,i+Ai + 1 - <i - <i+i) ( 7 - 2 ) 

multiplies an ordinary Hubbard hopping term ; only unprimed particles enter into the potential. 
There are four types of fermionic particles, respectively generated by cj i , ct i , i , c} i so that they 
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define a 16 dimensional vector space on each site. The corresponding numbers of particles are 
conserved. 

The factor M' aii+l works on a 4 x 4 space and its eigenvalues are ±1 with two- fold multiplicity. 
In particular this means that it cannot vanish, N' a i i+1 ^ 0. Moreover, if no primed particles are 
present, M' ai i+1 = 1 , V a, i. The same is true if the system is fully filled with primed particles in 
which case M' ai i+1 = — 1 therefore two of the sectors described by this Hamiltonian are equivalent 
to the ordinary Hubbard model. A Russian doll structure is appearing: if the projectors are well 
chosen, a larger model contains the small ones. 

If there are primed particles only, the energy vanishes (but not momentum). If the potential is 
interpreted as a Coulomb repulsion, then unprimed particles only carry electric charge. 

The compound objects formed by ci-c^^ are rigid: no other term in the Hamiltonian can 
destroy them. In this sense, we have four types of carriers, with the same charge but different 
behaviours: two are the elementary objects c' ai in two polarisations a =f , j, two are the compound 
objects (in two polarisations). 

The symmetry, according to is gl(l\l) © gl(l\l) © ^(1|1) © gKM 1 )- 

At second order in jj the following effective Hamiltonian appears 



2) 1 T 1 

i=i 



(7.3) 



that looks like a deformation of an XXX model. It has an enhancement of symmetry with respect 
to (|7.ip in the sense that its symmetry is gl(2\2) © gl{2\2). 

The two sites action of (17. 3p is a 64 x 64 matrix that can be easily disentangled leading to both 
the blocks given in (|5.5p . In summary, it has eigenvalues and 2, with multiplicity 48, 2 with 
multiplicity 16. 



7.2 gl(4) © gl(A) Hubbard Hamiltonian 

We consider the model based on gl(4:) © gl{4) and take the projectors according to the example 
(13. 3p with A?" = 4 , M = 0, in two copies (up and down) 

L 

H-Ruh = ^2 | ^2 { C t,i C a,i+l C },i C (T,i+l + C a,i+l C a,i C l,i+l C a,i + 

+ <,i<i+i(i,i c a,i+i + c l,i+i c *,i) + n«T,tn (T ,i + i(cJ ji t4 ji+1 + cj i+1 c^)) + 

+ U(l - 2n hl n' tl )(l - 2n Ll n' Ll ) } . (7.4) 

This model has the same vector space dimension of the gl(2\2) one (17. ip . dim(V) = 16 but the 
elementary projectors are different and lead to slightly different interactions. That the two Hamil- 
tonians are different is manifest if one examines the original form (|4.8|) with the basis matrices E a p , 
before the Jordan- Wigner transformation. 

Here there is complete symmetry between primed and non-primed particles; the effect of 
Coulomb repulsion only appears when both primed and unprimed particles are on the same site; if 
one of these types is alone, no Coulomb interaction is felt. Observe that \in! ai = \ everywhere and 
for all polarisations (or else if n a ^ = 1), we re-obtain the g/(l|l) Hubbard model. The kinetic term 
also has a strange feature: a particle is allowed to move only if it is accompanied by a particle of 
the same polarisation (i.e. up with up) but opposite type (i.e. primed with unprimed). 
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7.3 #/(4|4) © #/(4|4) Hamiltonian 

Following the example (|3.3p . the following Hubbard Hamiltonian is obtained 
L 

H Hub A) = Yl { Yl (4,i^,*+l + C U+l C °,i + 1 - ~ (4/*,i+l c *,i c ",i+l 
i=l <t=T,J. 

C o-,i+l C o-,i C o-,i+l C cr,i ~~ n a,i n (7,i+l\ C a,i C (T,i+l + C o-,i+l C cr,J 

" <i< i+ i(<i<i+i + + ^(1 " 2n' T ,n' T V)(l - 2n' ii4 n'/ i4 ) } (7.5) 

Here there are six types of fermions, c a , c& ,Co so the local (one site) space of states is 64 x 64. 

One observes that this Hamiltonian exhibits a 'Russian doll' structure. Indeed, there are four 
sectors in the space of states where the <?Z(4|4) Hamiltonian reduces to the gl{2\2) one, that also 
reduces to the </Z(l|l) one. For example, one sector is given by n'^ i = n\ i = 1 for 1 < i < L. 

8 Two-particles interaction 

We sketch here the preliminary effects that we observed studying two particles in interaction. It is 
convenient to consider a reference state as being a particle "vacuum" (pseudovacuum) 

ft = (el ® e[) ® (e{ ® e{) ® . . . (e| ® e{) (8.1) 

12 £ 

where index under the tensor product symbol labels the lattice sites. All other states are considered 
excitations of this pseudovacuum. Then particles are distinguished by the type, according to the 
subspaces W, W: 

a,b,... G W 

a, 6,... g w 

and an upper index |, J, will be added to distinguish polarisation^. 

Within the universal XX models, all particles satisfy the exclusion principle, namely they cannot 
appear on the same site. If two particles are both from W or both from W, they reflect each other; 
if they are one from W, one from W, they traverse each other by tunnel effect. 

In the universal Hubbard models, the coupling activates a sort of electrostatic interaction felt 
by particles of opposite polarisation only. Indeed, the potential term in (|4.8p squares to the identity 
(13. 6p so on one site states it has eigenvalues ±U. Which sign occurs is dictated by the membership 
to W or W according to the rule: with U > 0, equal type particles a)a^ or aJa^- repel each other 
but different type particles a^a^ or aJa^- attract each other. Observe that the vacuum itself is in 
the repulsive case so actually the only "visible" effect is the attractive one. 

9 Conclusions 

We have constructed universal XX and Hubbard model Hamiltonians based on general properties 
of projectors and permutations. The underling algebraic structure could be an ordinary or graded 

6 As already remarked, notice that particles of different polarisation or different type are not to be understood as 
conjugated: for example, a T and er are different objects 
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algebra gl(N\M) or possibly and infinite dimensional algebra. We have full control of the symme- 
tries of the models and we have performed the perturbative calculation d la Klein and Seitz [18] in 
the large coupling limit. 

We have emphasised that the gradation makes the Jordan-Wigner transformation a local iso- 
morphism. Therefore, the interpretation of the graded models in terms of 'electrons' is more 
natural. 

We discussed some examples, with their phenomenology. There, it would be very nice to see if 
the major screening effects observed (|7.ip and (|7.4p in the Hamiltonians can be interpreted in some 
condensed matter context. 

The next step in the study of our models is the determination of the spectrum and of the 
Bethe equations, as they were constructed for Hubbard or generalisation, using the algebraic Bethe 
ansatz [27-29,31] and the coordinate Bethe Ansatz of Lieb-Wu [11]. This is an heavy calculation 
which we postpone for further publication, but from the analytical Bethe ansatz approach, one can 
guess their form. In particular, as for spin chain models, one expects as many presentations of the 
Bethe equations as there are inequivalent Dynkin diagrams. All these presentations should lead 
to the same spectrum. For more informations, we refer to [32,33] where similar calculations were 
performed in the case of XXX super spin chains. 

Our models are graded by construction so they naturally contain bosonic as well fermionic 
degrees of freedom. We are working on examples with bosonic particles, that necessarily will be on 
infinite dimensional algebras. 

Finally, the Bethe equations will allow us to keep in touch with super-symmetric gauge theories, 
where integrability appears precisely in relation to the Hubbard model. 
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